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ABSTRACT 

Irrigation development is not an unmixed blessing. 

There are usually some undesirable environmental consequences 

such as waterlogging, salinity, increase in the incidence 

of malaria etc. It has been the experience in the past 

that there is a continuous rise in the ground water table 

when irrigation starts. An area is said to become water- 

-1 

logged when the water table is less than metres below 
ground level. Antiwaterlogging measures usually consist of 
drianage (surface and subsurface)^ embanl-cment and other flood 
protection works, installation of shallow tube wells and 
pumping of water to lower the water table, lining of canals, 

Y6 

suitable corpping pattern, land levelling^ etc. 

In the present work only one method of antiwaterlogging 
measure is studied, namely controlling of waterlogging by 
ground water management. A specific area namely hahi Right 
Bank Canal (MRBC) Command area is taken up for the study. 

While actual waterlogging on a large scale in this area 
has not yet taken place, its likelyhood exists, A mathematical 
model is developed for the aquifers of the FRBC Command area 
to study different policies of ground water management for 
antiv;aterl egging purpose. 

The mathematical model was validated using the 
readily available data. Sensitivity analysis on the 
parameters S,T,0 and L was made, It was found that S,T are 
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not very sensitive compared to Q and L, In order to vanity 
the effect of different boundary conditions on the results, 
three types of boundary conditions v^ere studied. It was 
found that the different boundary conditions studied do 
not have much affect on the results when we move away from 
the boundary. 

The validated model is then used for controlling the 
ground water table, One of the policy/^tested was to see 
how much pumpage is needed to lower the ground water table 
3 metres below ground level where necessary. It was found 
that about MCM of water per year is required to be pumped 
for this purpose. This waiter, of course, xs, available for 
summer irrigation, A second policy of lowering ground water 
table to 6 m below ground level was also studied. This 
yields about OSAO' MCM of waber per year. This additional 
water can be used not only in Summer but also for Rabi 
irrigation. 



CHAPTER I 


INTRODUCTION 


1.1 GENERAL 


Agriculture is a major sector in Indian economy 
accounting for almost half of the national income and provides 
employment to about two"-thirds of the work force. Accordingly, 
the greatest use of water resources has been for agricultural 
development which would benefit considerably by improvements 
of the management Water Resources' [Udai P. Singh, 1 984] . 

India’s agriculture is basically dependent upon the 
south-'west monsoon which is confined essentially to the 
period June to September, The rest of the year Is practically dry 
in almost all parts of the country. Therefore the development 
of irrigation facility is necessary through the conservation 
of the monsoon rainfall in surface reservoir and as ground 
water which should be utilized in the non-'monsoon seasons. 

There is enough literature to say that irrigation 
was practised extensively in ancient India. The accent on 
efficient water management in ancient India is seen in the 
following writings of Magasthenes, the Greek Ambassador in the 
Court of Emperor Chandra Gupta in about 300 B.C, [lARI Research 
Bulletion, 19B3 ]. 

' ♦The whole country was under irrigation and very 
prosperous because of two crops grown in a year with irrigation 
facilities. The District officers measure the land and 



2 


inspect the sluices by which water is distributed in the 
branch canals (water courses) so thaL every one may share 
his fair share of benefit''* 

It may be noted that irrigation plays a protective 
as well as productive role in the use of other agricultural 
inputs and thus becomes catalytic agent in enhancing the 
spread of ' 'green revolution' ' , 

1 .2 WATERLOGGING' PROBLEM 

Irrigation development, as v;ith any other development, 
is not an unmixed blessing* There at^e usually some undesirable 
environmental consequences such as waterlogging, salinity, 
increases in the incidence of malaria etc. It has been the 
experience in the past that there is a continuous rise in the 
ground water table when irrigation starts. For example the 
water table which was about 30 metres below the ground level 
before the start of irrigation in the Sindh area in the Indus 
River Valley (present Pakistan) started rising at a rate of 
about 1/3 metre per year and presently large areas are water-- 
logged and saline. The opening of Ganga canal resulted in 
rise of water ta le in the Ganga-Yamuna Doab from a depth 
of 12.2 metres to abouti4.6 metres in about 100 years* 

Similarly in the area commanded by the western Yamuna canal 
the wa’er table rose at an average rate of 16 cm per year 
during the years 1932 to 1963 (Irrigation Commission Report, 
1972). 
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An area is said to ‘become waterlogged when the water 
'] 

table is less than metres below G.L, High water table, 

heavy monsoon rainfall, perennial irrigation and flat nature 
of terrain all contribute to waterlogging. Anti waterlogging 
measures usually consist of drainage (surface and subsurface), 
embanlcment and other flood protection v'^orks (to check 
inundation of lands' , installation of shallow tube-wells, 
lining of canals (to reduce seepage), suitable cropping 
pattern (using crops needing light irrigation) , land levelling, etc. 
1.3 ANT I- WATERLOGGING MEASURE 

In the present work only one methoc? of antiwaterlogging 
measure iA studied, namely controling of waterlogging by 
ground water management. It is to be understood that in any 
given situation perhaps all the measures mentioned above will, 
have to be undertaken to a lesser or greater degree depending 
upon the situation, 

A specific area namely Mahi Right Bank Canal (MRBC) 

command area is tak.^n up for the study. While actual vmter- 

logging on a lar^e scale in this area has not yet taken place, 

its likelyhood exists. This area is selected because a lot 

•f published data is available for this area. The main source 

of 

of information is the publication /the V/ater Technology Centre, 
lARI: ' ’Resource Analysis and Plan for Efficient Water 
Management- A case study of the MRBC Command Area, Gujrat ,1 983 ’ * . 

Ground water was the major source of irrigation during 
the period prior to the commissioning of the MRBC irrigation 
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system. It continues to be so even after tvro decades of 
introduction of canal irrigation. The use of ground water 
has progressively increased over the years. In spite of 
this, the water table in the command area has been rising 
continuously at an ever increasing rate. V/aterlogging has 
become a recurring problem in the command area [lARI Research 
Bull. 1983]. 

In the present study a mathematical model is 
developed for the aquifers of the MRBC Command area to stv.''"y 
different policies of ground water management for anti- 
waterlogging purpose. 



CHAPTER II 


THE S'^UDY AREA 


2.1 LOCATION 

Mahi Right Bank Canal Command area is situated in 
the heartland of the proposed Narmada high level canal 
command area in Gujrat which has an estimated potential of 
over two million hactares (see Fig, 2,1), 

Kaira district, in which the MRBC command area is 
located, is one of the 19 districts of Gujrat State, It 
lies between 22°?^ and 23^13' north latitudes and 72^25’ and 
73^37* east longitudes* The geographical area of the district 
is 6888 sq, km,, which is about 3,67 percent of the total 
area of the State. 

2.2 GEOLOGY 

According to geologists, the MRBC Command area is 
covered by the alluvial deposits carried by Mahi and Sabarmati 
rivers. The area is the result of the continuous degradation 
of the main rocks namely, phyllites, schists and quartzite. 

The alluvial deposits overly the Deccan trap basements. 
The total thickness of alluvial deposits which are unconsolidated 
varies between 30 to 150 metres from east to south-west, 

A typical cross section of the area is shown in the 
Figure 2.5. It indicates that the aquifer in the MRBC command 
anea consis-te of gravel and sand with an average thickness of 
about 30 metres and varying upto 60 metres, A distinct 
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water table aquifer of about 10 tu 1b metres thickness over- 
lying the main confined aquifer separated by layers of clay 
and kankar is seen. In general clays are mixed with kankar 
and sandy material. 

2 . 3 HYDRO-METEOROLOGICAL CHARACTERISTICS 

2.3.1 RAINFALL 

Indian agriculture is largely dependent on monsoons 
and it greatly suffers due to erratic distribution of rainfall 
over the year and its unpredictable behaviour. The faf^ure 
of monsoon in a year or its late arrival considerably affects 
the economy of the nation [Singh, R.L. 1978], 

The average annual rainfall of the MRBC Command area 
is 823 mm, spread over 35 rainy days . 96 percent of the 

annual rainfall is received during south-west monsoon. 

The Indian Meteorolo^iical Department maintains the 
records of rainfall, A record of monthly rainfall for the 

year 1976-79 is given in Table 2.1 . 

/)/ 

2.3.2 DRIANAGE 

The Mahi Right Bank Canal Command area is characterised 

by a generally flat topography, restricted natural drianage, 

comparatively heavy soils and semi-arid climate, all of vrhich 

of 

are conducive to the development^water logging and salt 
accumulation [lARI Research Bull. 1983]« 

The possible causes of the continuing problem of 
water logging are attributed to the following; 



TABLE 2.1 : MONTHLY RAIl^FALL IW THE STUDY AREA DURING THE YEARS 1976-1979 
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1 . Inadequate number and capacity" of drains to effectively 

dispose off the excess rain water from the cropped 
fields. 

2, Improper irrigation in the command area. 

3 • Seepage from canals system and consequent rise of 

water table . 

4. Inadequate maintenance of the drains. 

A systematic approach on the provision of drainage 
facilities in the MRBC command arfra was initiated in 1971. 

The details of the drianage systejn are given in Fig* 2.2. 

2,3.3 CLIMATE 

2 . 3 . 3 . 1 TEMPERATURE AND HUMIDITr 

The climate of Mahi Righ"^ Bank Canal Command area is 
semi -a rid. It is characterised by hot summer and general 
dryness, except during the south west monsoon season which 
experiences heavy rain. It is divided into four distinct 
seasons, cold season from Dec, to Feb, is followed by hot 
summer from March to middle of June* Mid June to Sept, is 
south-west monsoon season, October and November post monsson 
season, 

India Meteorological department has not maintained 
meterological observatory within the command area for recording 
temperature, humidity, sunshine and other vi^eather data. The 
information on temperature presented here is based mainly on 
the records of the observatory at Baroda close to MEBC area. 




Fia2'2 Main drains provided in the study area 
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During the period from March to May there is continuous 
increases in temperature. May is generally the hottest month, 
the mean daily maximum temperature being 41 and mean daily 
minimum temperature about 26®C, There is appreciable drop 
with the onset of monsoon season* January is generally the 
coldest month with the mean daily maximum temperature of 
about 30^C and mean daily minimum temperature about 11 
During the cold waves, the temperature sometimes drops to 
2 to 

Humidity is generally high during the south-west 
monsoon season, April -to June is the driest part of the 
season when the relative humidity is 46 to 91 percent in the 
morning and 1 6 to 75 percent in the afternoon. 

2. 3. 3. 2 WIND 

During most of the monbhs in the year, the winds are 
generally light, ranging from 2.6 knv'hr to 9.8 kny^hr. 
Comparatively strong wind blow in summer months, 

2.3.4 PAN-EVAPORATION 

Pan evaporation values for the study area are not 
available. However, some values are available for a neighbouring 
area, namely the Vasad area. Table 2.2 gives the pan evapora- 
tion data, observed at Central Institute of Soil and Water 
Conservation, Vasad [lARI Research Bull, 1983 ]. 



Table 2.2 ; data of average daily evaporation of vasad 
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2.4 WATER TABLE 

The ground water levels in the MIBC comrnand area are 
being observed since the introduction of canal irrigation. 

The water level data of open wells and shallow tube wells 
were collected twice in a year to represent the pre-monsoon 
and post-monsoon water levels. Based on the analysis of pre- 
and post-monsoon water levels for the period 1961 to 1973, 
it was observed that on an average^ the water level rise was 
0.28 metres per year [lARI Research Bull. 1935]. There are 
some regions in the command area where the water table rise 
was observed to be much steeper. The trend of water table rise 
at selected locations is shown in Fig, 2,3 [lARI Research 
Bull. 1983], 

It may be seen that by 1980 the water table rose 
considerably in the northern and north western parts of the 
command area as compared to the pre-irrigation period. The 
changes in the depth to water table in the same region during 
the pre and post monsoon period in the year of 1980 are 

f 

presented in Fig, 2,4. 

Table 2,3 gives the area with different ranges in 
depth to water table during pre-monsoon periods in 1958 and 
1981 in the MRBC command area (upstream of Alang drain) and 
Table 2.4 gives the areas under different groups of depths 
to water table during pre-monsoon and post-monsoon periods 
for the year 1980, It is seen from these that while the 
waterlogging may not be acujse at present (1980), it is on the 




Fig.2'3 Depth to water table (premonsoon) 




Fig. 2-4 Pre and post monsoon wat 
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TABLE AREAS WITH DIFFERENT RANGES IN DEPTH TO WATER 

TABLE DURING PRE-MONSOON PERIODS IN 1958 At>ID 1981 


IN 

THE MRBC 

COMMAND AREA 

(UPSTREAM OF ALAMO DRAIN) , 

Depth to water 
table from 
ground surface 


Area 

covered 



1953 

1981 


m 

ha 

Percent 

ha 

Percent 



of total 


of total 

1 .5 

- 

- 

475 

0.20 

1,5 to 3.0 

2,643 

0,90 

9,360 

3.20 

3.0 to 6,0 

18,053 

6.10 

83,333 

28.30 

6.0 to 9.0 

1 5,277 

5.20 

33,500 

11,40 

9.0 

2,57,382 

87.80 

1,67,207 

56.90 

Total 

2,93, 860 

100.00 

293,860 

100,00 


Source : Report on water conditions in MRBC Area, Soil Survey 

Division, Dept, of Irrigation, Govt, of Gujarat, 1 931 . 
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TABLE 2,4 : AREA UI^ER DIFFERENT GROUPS OF DEPTH TO WATER 
TABLE DURING PRE-MONSOON AND POST-MONSOON 
PERIODS FOR THE YEAR 1980 IN THE STUDY 
COmkm AREA (UPSTREAM OF ALANG DRAIN) 


Depth to water 


Area 

covered 


table from ground p 

surface 

-monsoon, 1 980 

Post-monsoon, 1 980 

m 

ha 

Percent 
of total 

ha 

Percent 
of total 

1 .5 

250 

0.10 

10,725 

3.70 

1,5 to 3,0 

7,780 

2,60 

'll ,335 

1 7.50 

3.0 to 6,0 

71 ,025 

24.20 

48,911 

1 6 . 60 

6.0 to 9.0 

hZ , 077 

U.30 

29,009 

9.90 

9.0 

1 ,72,728 

58. SO 

1 ,53,830 

52.30 

Total 

2,93,860 

100.00 

2,93,860 

100.00 


Source: Report on water table conditions in MRBC Area. Soil 

Survey Division, Dept, of Irrigation, Govt, of 
Gujarat, 1981, 
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TABLE 2.5 : MONTHLY RELEASES OF IRRIGATION WATER AT THE 

WANAKBORI WEIR HEADWORKS MRBC SYSTEM (BASED ON 
DAILY OBSERVED VALUES). 




SI. 

Month 

1 976-77 

1 977-78 

1978-79 1979-80 

1980-81 

Mo. 


ha-m 

ha-m 

ha-m 

ha- m 

ha-m 

1, 

June 

1 ,410 

0 

0 

0 

0 

2. 

July 

16,972 

9,374 

16,846 

1 6,662 

24,747 

3. 

August 

8,944 

24,433 

31 ,683 

26,331 

29,236 

h. 

September13,933 

26,995 

33,387 

33,202 

45.667 

5. 

Oct. 

23,295 

26, 640 

32,352 

25,371 

35,430 

6. 

Nov, 

8,573 

11,698 

10,883 

6,475 

9,609 

7. 

Dec . 

6,723 

9,996 

4,362 

3,143 

8,623 

8 . 

Jan, 

2,224 

7,680 

6,894 

5,751 

11 ,166 

9. 

Feb, 

1,798 

4,738 

7,617 

5,584 

11 ,126 

10. 

March 

1 ,055 

2,70'5 

3,019 

5,577 

11 ,199 

11. 

April 

0 

0 

0 

7,005 

0 

12. 

May 

0 

0 

0 

4,560 

0 

Total 

Year 

84,927 

1 ,24,189 

1,47,043 1 

, 44 , 661 

1 ,84,847 

Average annual 

flow for 

the period 1976-77 to 

1979-80: 

1252 MSM 

Source: 

Department of Irrigation, 

MR3C, CAD A, 

Nadiad, 
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increase as more and more area is coming within the water 
table depth of about 0 to 1,5 metres below ground level. 

Hence suitable measures require to be taken to control this. 

One of the metnods is proper ground water management, 

2.5 IRRIGATION 

Irrigation water is draivn from the MRBC irrigation 
distribution network as well as from the ground water basin. 
Irrigation from other independent sources, such as tanks, 
local streams etc, is negligibly small. The total volume of 
irrigation water delivered through the MRBC distribution 
network during period June to May for the years 1976-77 to 
1979-80 are dhown in Table 2,5, 

The area irrigated by wells during the years 1976-77 
to 1979*“60 for the period June to May is presented in Table 2,6, 

TABLE 2.6 : AREA IRRIGATED AND DRAFT FROM WELLS IN MRBC COMMAND 
AREA 


Year 

Kharif 

ha. 

Draft 

mm"^ 

Annual area 
ha 

Draft 

1 976-77 

53411 

67 

91 794 

104 

1977-78 

54075 

68 

1 07002 

121 

1 978-79 

66822 

m 

133510 

151 

1 979-80 

69371 

87 

120575 

136 

Average 

60920 

76 

113220 

12B 


+ Depth of water distributed over the GCA of 31600 ha 
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It is evident that at the ’oresent rate of development 
in irrigation activity, the ground water table will continue 
to rise and will result in waterlogging. It is therefore 
necessary to control it. Management of groundwater through 
proper development of ground water will change the ground 
water regime and the incipient conditions of waterlogging 
are controlled. The exact Impact of ground water development 
in different regions can be assesed by formulating a computer 
based ground water management model. 

2.6 AQUIFER CHARACTERISTICS 


In the command area pump-test analysis has been done 

in several wells, these tests shows that coefficient of 

transmissivity lies between 186 m /day to 6830 m /day with 

the higher values occuring in the eastern parts and lower 

values in the western parts. Long duration pumping tests 

conducted at selected sites Indicated that the storage 

coefficient values of the acjuifers In the MRBC command area 

-A “3 

range between 1.5x10 to 4.17x10 . These tests also .indicated 

that the aquifer is leaky in most of the area. The average 
values of leakage coefficient (k’/b’ ) was observed to be 
'^tOOI' 6' day , Table 2,7 gives depth, thickness and coefficient 
of permeability of aquifers at various locations. 
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TABLE 2 . 7 : 

RANGE OF GEOHTOROLOGI" 

TALUKS COMPRISING THE 

\L PARAMETERS 

MRBC COMMAND 

IN THE DIFFERENT 

AREA 


Taluk 

Depth 

Thickness of 

Hydraulic 




aquifer 

conductivity 



(m) 

(m) 

(m/day) 

1 . 

Bor sad 

91 -1 67 

26-62 

40-105 

2. 

Thasra 

40-65 

10-22 

80-230 

3. 

Anand 

52-1 92 

13-48 

34-88 

4. 

Cambay 

1 1 0-1 54 

21-60 

31-70 

'5. 

Nadied 

134-152 

27-61 

30-75 

6 . 

Petlad 

97-1 66 

28-70 

3S-85 

7. 

Matar 

1 52-1 55 

18-40 

21-52 


CHAPTER III 


GROUND V/ATER MANAGEMENT MODEL 

Several methods are used by ground water hydrologists 
to tackle ground' water problems. Some of these are (1) Analytical 
methods (2) Analog methods (3'^ Digital methods and (A) Sand 
model methods etc, 

3 .1 ANALYTICAL METHODS 

Analytical methods are used for aquifers having 
uniform properties. These properties, in general, are 

1 . Coefficient of transmissibility 

2. Coefficient of storage 

3 . Coefficient of leakage 

The mathematical equations used for simplified 

aquifer systems are given in Table 3.1. 

TABLE 3.11 SIMPLE ANALYTICAL MODEL EQUATIONS 


SI, No. Mathematical Model Equations Applied to 


2 , 


Q 


s 


U 


A^T 

r^S 


W(u) 


U “ 


ATt 

Q 

r^s 

Trt 


w(u —-,y) 
' m 


Isotropic nonleaky 
artesian aquifer with 
fully penetrating wells 
and constant-discharge 
conditions 


Isotropic nonleaky 
artesian aquifer 
with partially 
penetrating wells and 
constant discharge 
conditions 


Contd 
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Table 7 contd,.. 


Si. No* Mathematical Model Equations 
used in metric system 


3* 



m 



4T t 


Applied to 


Isotropic leaky 
artesian aquifer with 
fully penetrating wells 
and constant discharge 
conditions without water 
released from storage 
in aquitard. 




^ Q 


Ko (f') 


Q . r 

^ ^t 


u„ * 
a 


u 


r^S 
4 T t 

4 T t 



Isotropic water- 
table aquifer with 
fully penetrating wells 
and constant discharge 
conditions . 


where, 

3 


dradown in meters, 
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Q « discharge in m^/day 

T » coefficient of transmissibility of aquifer 
in m^/day 

S - coefficient of storage of aquifer, fraction 
r = distance from production //ell to observation 
point in metres 

t = time after pumping started in days 
m = saturated thickness of aquifer in meters 

= distance from top of aquifer to top of screen 
in meter, 

K ‘ = coefficient of permeability of aquitard in 
m/day 

= saturated thickness of aquitard in meters 
Sy = specific yield of aquifer tn meters 
Kq 5= Besselfunction of the zeroth order and second 
kind, 

W(u), V/(u,r/B), W(u,r/Dt), W(u,r/B,i^) are well functions* 

Principle of superposition is used in applying the 
analytical methods mentioned above. Since the approximate 
governing differential equation used in the simplified models 
is linear, the effects of pumping of or recharge from several 
wells are superposed, Tf the number of wells operating in an 
aquifer is small, this addition can be done manually. If th^ 
number of wells is large then recourse is taken to digital 
computers , 
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Boundary conditions are usually taken into account 
in these analytical methods by the method f images [Walton, 

1 970] . 

3.2 DIGITAL MODEL-LITERATURE REVIEW 

Aquifers encountered in nature are generally complex. 
The aquifer properties vary from place to place. Such aquifers 
are called non-homogeneous. Further, some properties such as 
aquifer transmissibj lity may vary depending on the direction. 
Such aquifers are called anisotropic. The boundaries are 
generally irregular and the boundary conditions are not 
generally well defined. It is very difficult to apply simple 
analytical methods in such situations. Numerical methods 
are generally used along with digital computer to solve such 
problems . 

With the widespread availability of digital computers 
has come the development of mathematical models of aquifers. 
Applications are expanding, programming techniques are steadily 
improving, and computer capabilities are growing so that it 
is safe to say that almost any type of ground water situation 
can be studied by means of a digital computer model [Kleinecke, 

1971] . 

There is voluminous literature on mathematical 
modelling, analog and digital computers applied to ground 
water problems. 
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Prickett and Lonnquist (1968) gave a comparative 
study of analog and digital simulation technique applied 
to various situation, 

Pinder and Bredehoeft (1968) discussed the application 
of the digital computer for aquifer evaluation » 

Rushton and Bannister (1970) reported on slow time 
resistance capacitance analog as an alternative to usual 
fast time analog. 

Lakshminarayana (1971) reviewed electric analog 
model for management of aquifers with the help of some 
examples of mathematical models for simple and complex 
geometric boundaries, According to him the models can be 
useful if the historical data are reliable. Reliable analy- 
tical, analog and digital models can then be developed for 
efficient management of the aquifers. 

Cooley and Sinclair (l 976) investigated some of the 
methods and conditions necessary to determine bounds on 
values of the hydrogeologic parameters which characterise 
a model of steady state area near surface ground water 
flow. 

Richard L, Cooley (1977) developed & nonlinear least 
squares solution for the hydrogeologic parameters, sources 
and sinks, and boundary fluxes contained in the equations 
approximately governing two-dimensional or radial steady 
state ground water motion, through use of a linearization 
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and iteration procedure applied to the finite element 
discretization of the problem, 

Trescott, P,C, and S,P, Larson (1977) compared the 
efficiency of line successive over relaxation with a two- 
dimensional correction procedure (2DC), the iterative, 
alterating direction implicit procedure (ADI), and the 
strongly implicit procedure (SIP) to solve finite difference 
equation used to simulate several ground water reservoirs. 

McElwee, D. Carl and M. Arif Yukler (1978) studied 
the systems response to various parameters. It is shown here 
how to evaluate the perturbed hydraulic head for a small 
change in aquifer parameter by means of a first-order 
sensitivity formation, 

Prickett (1979) addressed in his paper both the pros 
and cons ground water modelling and presented from a 
neutralists stand point, 

Cooley, R,L, (1979) applied techniques of nonlinear 
regression to estimate the hydrologic parameters (values of 
hydraulic conductivity or transmissibity) , recharge, discharge, 
and boundary fluxes for steady ground water flow models of 
two field areas; Truckee Meadows, Nevada, and a cross section 
in the Hula Basin Israel, 

James and Charles (1980) gave an over view of ground 
water modelling. The same authors discussed the different 
approaches to be used in the model. They said that for any 
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given class of problems the choice of the best approach 
depends on the processes being modeled, the accuracy desired 
and the effort that can be expanded on obtaining a solution. 

They in the paper represented the use of a ground 
water flow model to analyse an aouifer system comoosed of 
glaciofluvial deposits » 

Munter and Anderson (1 981 ) developed a model to 
estimate the additional data ’'seepa'^e rate’’ foi‘ the 
Has Lake and Nepco Lake, Wisconsin. 

Kant, Naney and Witz (1982) demonstrated the 
application of a predictive ground water potentiometric 
head model to estimate the profitability of irrigation in 
constrast to that of dry land forming, 

Keating (1982) used a lumped parameter model to 
explain an usual chalk aquifer consisting of a thin zone 
of high storage and transmissivity whose response of small 
annual fluctuations and high stream flows were considered 
to be anamolous. 

Chu, Wennsen and Robert (198 h^) px'toented a simpler 

alternative which permits an easy direct solution of the 
Bousnesque equation. They used a forward in time, central 
in space (FTCS) explicit finite difference method in the 
simulation model. 
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Paul Richard Eyre (1935) simulated the flow of ground 
water and the effects of future water development using 
AQIFEM, a two dimensional finite element flow model, modified 
for a^luifers containing a sea water interface in south- 
eastern Oahu, Hawaii. 

3.3 MATHEMATICAL MODEL 


In applying a distal model the aquifer is first 
descretized by superimposing a grid of horizontal and vertical 
lines in the finite difference method^* The intersection 
of the grid lines is termed as nodes. 


The flow equation is then used for each node. The 
flow equation is based on the equation of continuity and 
Darcy’s law. The equation of continuity can be written as. 
Inflow - Outflow « Change in storage . 

This relation with Darcy’s law yields the equation 


where , 








r 


T » transmissibility of the aquifer in the x and y 
direction in the horizontal plane, 

S s storage coefficient 

q = net ground water withdrawal per unit area 
r * leakage into the aquifer per unit area 
r 

L ^ rr “ h) ] 

b’ 



33 


K’ « coefficient of permeability of the overlying 
semipervious aquifer (m/day) 
b* a the thickness of the overlying semipervious 
aquifer (metre) 

H a head in the overlying aquifer (metre) 
h = head in the main aquifer 

X (5< y = rectangular coordinates 
t « time coordinate* 

For each node this equation is written. The 
differentials and "^y and 7)t are approximated by 

differences Ax, Ay and At - The discretized model is a 
reasonable representation of the contiuous system provided 
Ay are small compared to total area of the aquifer. 

The difference equation for each node can be derived 
either from the differential equation given above or from 
physical considerations using Darcy ^s law and continuity 
equation. In the implicit method the space derivatives are 
replaced by finite differences at the time at which the heads 
are to be calculated. This involves a set of simultaneous 
equations but the advantage in this method is that it is 
unconditionally stable regardless the size of the time 
increment. 

Using an implicit finite difference scheme for the 
present problem, we get the following Equations (Prickett,1 97--^; . 
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( ) + T. . CL, ) 

' > U ^ ._2 o . o 




AJ^ 


j . -1 


<Ax 

h^- .-h^ 
At 


Q 


R. 


ilj -S , , t\ . . 

s S. . ^ JiAI *. - n 

Ax Ay Ax^y 


where, 


AX, Ay 


i.O 

At 





id 

i.d 

i.d 


finite difference grid size in the x and 

y directions (metre) 

column and row index of a node (i,d) 

time increment elapsed since last calculation 

of heads(day) 

calculated head at the previous time 
increment at node (i,d) (metre) 
calculated head at the present time 
increment at the node (i,o) (metre) 
transmissibility of the branches between 
nodes (i,j) and (i**-!,^) and between nodes 
(t,j) and (i,j+l). The transmissibility of 
a branch between any two nodes is taken to 
be the average of the values at these two 

p . 

nodes (m /day) 

storage coefficient , at the node (i,j) 
(dimensionless) 

net withdrawal rate from the element of the 
aquifer centered at (i,j)(m^/ day) 
leakage into the element of an aquifer cantred 
at (i , j)(m^/day) . 



- K ' 

.flF /— = AxZ\y = Leakage factor, 

b ' 

The equation is written for each node which results in a 
set of simultaneous equations. 

Many methods are available for solving a set of 
simultaneous equations. Some of these methods are the 
successive over relaxation method, alternating direction 
implicit method, and others. The simplest one is successive 
over relaxation method, hut it is adequate for small size 
problems. In the present problem a modified version of 
iterative alternating direction implicit method (Pricket, 1 975) 
is used. 

A computer programme written by Mr. T.A, Prickett 
and C,G. Lonnquist of Illinois Stat^^'^ater Survey, Urbana, 
Illinois was adapted for the problems encountered in the 
MRBC Command area. 



CHAPTER IV 


ANALYSIS AND DISCUSSION OF RESULTS 

4.1 INTRODUCTION 

The work done for the control of waterlogging for 
Mahi Right Bank Canal Command area is presented in this 
chapter. Table 4.1 presents the history matching to obtain ^ 
the estimates of S,T,L and Q parameter value5for particular 
element represented by a node. In order to make sensitivity 
analysis the basin is divided into three zones. Table 4.2 
shows the results for the average head values for average 
values of the parameters. In order to Understand the effect 
of different boundary conditions, the sensitivity analysis 
of the parameter is repeated for different boundary conditions. 
Table 4.9 presents the situation of water table below ground 
level and pumpage required to bring down the water table 
below ground level by 3 metres and 6 metres. All the work 
done is discussed below in detail, 

4.2 HISTORY MATCHING; 

The total area of the Mahi Right Banlc Canal Command 
area is 5*6x10*^ sq.km. It is covered by a network of 13x11 
grid lines(Fig, 4.1). Grid spacing taken here is 6000 m. At 

n 

each node a well is assumed covering 3.6x10 sq. m. area. For 
matching historical data several years of data are required for 
the reliability of the models. Unfortunately only two years 
data were available , That is the heads (with respect to mean 
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seal level) of 1979 and 1980 were available. Approximate 
values of transmissivity, storage coefficient, leakage etc, 
were also available. Based on this information, the 
parameters, mentioned above, were roughly fixed for each 
node. Taking 1979 heads as input, the heads after one year 
(for 1980) were computed with the help of computer. It is 
obvious that in first shot the 1980 heads value (computed) will 
not be matching \«/ith the 1980 heads value (observed). So, for 
matching the history data the values of parameters (S,T,Q 6< L) 
were varied in a fixed range. The method adopted here is 
just trial and error and several runs were taken for matching 
the computed 1980 heads with observed 1980 heads. Table A.1 
represents the heads of 1979 (observed), 1980 (observed) 

and 1 900 computed for each node , 

In order to study the results to changes in the input 
parameters, senjtivity analysis is required. The Mahi Right Bank 
canal command area being large it was decided to divide this 
into 3 zones based on S and T values and mal<e the sensitivity 
analysis for these zones. 

Table A. 2 represents the average head observed (1979 aud 
1980) and head computed. These heads are average values worked 
out zone^ise from Table 4.1 . Using the mathematical model 
zone wise heads were also recomputed by inputting into the model 
zonewise average parameters. These are shown in Table 4.3. 

It may be seen that there is slight difference in the computed 
heads for the year 1980. 
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4.3 sensitivity analysis 

For each zone the average valuo of S,T,Q and L is 
computed. Taking these values and average head 1979 (observed) 
as input the heads are computed for 1980. These computed 
heads are compared with the observed 1980 heads. RUN-A 
represents computed heads for 1980 for the average conditions 
of parameters. This is used for the purpose of making 
comparision with other runs. The procedure adopted for 
sensitivity analysis is to change one parameter at a time 
holding all other parameters at their average values. The 
changes are ma‘.e within a reasonable range. RUN 1 and RUN 2 
are taken after decreasing or increasing the transmissivity 
values hy 50 percent. In RUN 5 and RUN 4 storage coefficiont S 
is decreased or increased by one order of magnitude of the 
average values. Similarly in RUN 5 and RUN 6 0 (pumpage ) is 
decreased or increased by 50 percent and in RUN 7 and RUN 8, L 
(leakage) values are decreased or increased by 50 percent of 
their average values. 

The boundary conditions imposed on the model are 
constant head boundary conditions for evaluating the average 
parameter values. It is felt that these boundary conditions 
may not be appropriate all along the boundary. At some places 
along the boundary it is possible a different type of boundary 
condition than a constant head boundary condition may be more 
appropriate. Hence it was decided to study the effect of 
different types^boundary conditions on the heads inside the 
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boundary. Three types of boundary conditions were 
investigated, 

1 . Constant Head Boundary: - The head is constant along the 
boundary. This condition is easily imposed by making the 
storage coefficient S assume a very large value on the 
boundary, 

2, No Flow- “ This condition assumes that no flo\^ takes 
place across the boundary. This condition is fulfilled by 
taking the transmissivity (T) values as zero on the boundry. 

3. Specified Outflow:- This condition assumes that a specif led 
outflow takes place across the boundary^ A part to the 
boundary is assumed to be impervious and on the remeining part 
is allowed to take place across it. This condition is 
easily imposed by assigning the quantity that flows out across 
a segment of the boundary as a pumpage assigned to the node 
nearest to the segment, and the inflow across the segment of 
the boundary as recharge assigned to the node nearest to that 
segment. 

Sensitivity analysis can also be useful for site - 
specific applications to indicate what additional data is to 
be determined and area where additional data are needed, 

4.4 CONTROL OF WATERLOGGING 

Waterlogging occurs when the water table is zero to 

1.5 metres below the ground level* Prom the model it is 
possible to compute^ the year 19B0 head© using the 1979 heads, 




These heads are measured with respect to mean sea level. 

We "know the geohydrological conditions of the aquifer at 
different section taken through the MRBC Command area. The 
ground levels are alsd known with respect to MSL, Computed 
heads can therefore be used to find the depth of water table 
below ground level. In the Table 4.9 the water table below 
the ground level is shown at each node. It is seen that 
at many nodes waterlogging has occurred. It is necessary 
to keep the water table below the ground level at a reasonable 
depth. One policy can be to keep the water table, say, 

3 metres below the ground level. The additional pumpage 
required to be done in this case is shown in the Table 4.9* 

The total ground water required to be pumped is about 
544 MCM per year. This extra water is available for 
irrigation, 

A second policy can be, from the management point of 
view, to stabilize the water table below ground level at 6 
metres. This will yield an additional quantity of about 
MCM per year which can be used for irrigation, 

4.5 DISCUSSION OF RESULTS 
4.5.1 MATCHING OF HISTORICAL DATA 

The validation of the model is done by matching 
historical data. For proper validation of ground water models 
at least S to 1 0 years of data are required, Th^n half the 
number of years of data is used for matching and the other 
half for checking. In the present case only two years data 



42 


were available hence the results can not be fully relied 
upon. 

While matching with historical data it is assumed 
that the heads on the boundaries are constant. But the 
available data show that appreciable changes In heads occur 
on the boundaries . For this reason the boundary nodes are 
omitted in sensitivity analysis. The results of history 
matching are shown in Table 4.1 node by node. It is seen 
from the results that inside the boundary the heads of 1980 
Oobserved) are fairly close to the 1980 heads (computed). 

The var’ious parameters for each of the nodes are shown in 
column (3 ), (4), (5) and (6). 

As mentioned before for making sensitivity analysis 
on the parameters the area in divided into three zones. For 
each zone average values of S, T, Q and L are computed and 
shown in Table 4.“^. Further using the mathematical model 
th? heads computed for each 35one for the zonal average 
parameters are also shown in Table 4.5. stated above 
the difference between these heads is small in each zone. 
These average heads will be used for purposes of comparision 
later in sensitivity analysis. Also sensitivity analysis 
has been carried out for different boundary condjtions. 

4.5.2 SENSITIVITY ANALYSIS 

Table 4.4 shows the results of sensitivity analysis 
for constant head boundary conditions. In this Table RUN 1 



represants results for T ^ 0t5 ,and all other parameters 
a»"e held at their average values. Similarly for the other 
RUNS the changed value of parameter is shown under the 
particular RUN in the table. 

By studying the results shown in this table (Table 4,4) 
the following cone lusiore may be drawn, 

(i) Effect of Variation in T:-- The results indicate that 
zone I and zone II are not affected with variation in T. 

Zone III is only slightly affected, 

(ii) Effect of Variation in S:- The results got after varying 
the storage coefficient by one o’-'der of magnitude above and 
below the average values show that there is not much change 

in Zones I and II whereas there is only a slight change of 
heads in zone III, 

(iii) Effect of Variations in Q:- All the three zones 
are sensitive to changes in U. Zone II is more sensitive 
and Zone I is least sensitive. 

(iv) Effect of Variation in L: Zone II is more sensitive 
than Zone III, Zone I is the least sensitive. 

Thus it may be stated that in general the parameters 
Q and L are more sensitive than the parameters S and T, 

Further Zone II is more sensitive to changes in Q and E 
comparedto Zone III and Zone X is the least sensitive. Thus 
more care and effort should be exercise in collecting data 
on Q and R than on S and T, 



kk 

(v) Effect of Boundai'y Conditionst- In order to see the 

sensitivity of boundary conditions, different boundary 

conditions were considered. The results for the three 

different boundary conditions are shovm in the Tables 4*5, 

5* 7 

4.@, 4.0. It is seen from the results that changes in 
boundary conditions do not affect the compu'ced heads inside 

' ' ■ ' MM , 

the boundary much. This is because the grid size is fairly 
large (6000 m x 6000 m) because of which the inside points 
are far removed from the boundary and hence changes in 
boundary conditions are not felt at these points. This would 
therefore mean that one can perhaps impose very simple 
boundary conditions such as no flow or constant head boundary 
conditions (instead of more complex boundary conditions such 
as mixed flow condition ) and still get acceptable results 
as long as the boundaries are far away from the point of 
interest. 

4.6 ANTIWATERLOGGING EFFECT 

As has already been mentioned one of the antiwaterlc'r' 
ing measure is by proper management of ground water. Having 
developed a mathematical model for the study area including 
validation of model and sensitivity analysis, we can say that 
we have an operational model for the purpose of ground v;ater 
management to centrol waterlogging. Using the model we can 
determine how much pumping should be done at each node to 
see that water table is lower down to say 3 m below ground 
level. Table 4.9 shows these results. It is seen from this 
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Table that the net quantity of water that is required to 
be pumped from the entire area comes to about 454 
This can be considered as extra ground water potential 
that is available for irrigation. As an alternative 
policy, we may decide that instead of lowering the water 
table only to 3 m below ground level we may lower it to 
6 m below ground level. This will not only remove water- 
logging but also give much more additional ground water. 
From the table it is seen that for this policy the net 
water available is about 3960 MCM per year. It is 
emphasized that this net quantity of pumpage should be 
over and above the pumpage required to account for 
(i) the fraction of applied ground water returning back 
as recharge after irrigation and (ii) the reduction in 
non-beneficial evapotranspiration due to lower ground 
water table. 


s 



TABLE 4.1 : RESULTS OF MTCHINS HISTORICAL DATA 
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TABLE 4,2: ZONSWISE 

AVERAGE VALUES 




Heads 1979 

Heads 1 980 

Heads 1980 



(Observed'^ 

(Observed) 

(Computed) 

Zone 

I 

1H.11 

18.68 

18.59 

Zone 

II 

9.59 

10.99 

10.01 

Zone 

III 

29.'=^? 

30.75 

30.83 


I 



TABL'S 4.3 : RUN- 4. MODEL RESULTS FOR AVERAGE VALUES OF FARAFETSRS FOR CONSTANT HEAD 
BOUNDARY CONDITION 
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TABLE 4.4: SEMSITIVITY ANALYSIS FOR CONSTANT HEAD BOUNDARY CONDITION 
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TABLE L.6: SENSITIVITY ANALYSIS FOR NO FLOV/ 'SOUNDARY CONDITION 
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^ Case 1 - V/ater table depth to be kept 3^ below ground level 
Case 2 - Water table depth to be kept 6 m below ground leve! 
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CONCLUSION AND SUGGESTIONS FOR FURTHER V/ORK 
•5.1 CONCLUSION 

Irrigation when practised improperly leads to 
several undesirable results such as waterloggingy salinity ^ 
incidence of maleria etc. Present thesis deals with the 
aspect of waterlogging. Several anti waterlogging measures 
such as drainage (surface and subsurface), embankment and 
other flood protection works (to check inundation of the 
lands), installation of shallow tube wells, lining of 
canals (to reduce seepage), suitable cropping pattern (such 
as switching over to crops which need light irrigation in 
preference to* crops which need heavy irrigation), proper 
land levelling (to avoid accumulation of water over land), 
etc,, are available. Many or all of these measures may 
have to be used in any given situation. The present work 
concentrates only on one method namely control of waterlogging 
by installaiion of tube wells and management of ground water. 
A mathematical model is formulated for this purpose. This 
model is applied to a specific study area namely Mahi Right 
Bank Canal Command Area, 

Published data available for the MRBC Command area 
was analysed. In this area after introduction of canal 
irrigation about 30 years ago there has been, a progressive 
increase in the ground water elevation. Many areas are 
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gradually coining Into the range cT waterlogged situation 
(when the water table below ground level comes to 0 to 1 
metres below ground level, the area is said to be waterlogged). 
Efforts have been made by the project authorities to construct 
drainage network, The progressive increase in the water- 
lofeg^d area indicates that this is not enough. Along with 
drainage other methods such as installation of tube wells 
and large scale pumping of ground water are also required. 

The mathematical model developed for the study area 
was validated using two years of ground water data readily 
available. It may be noted that one neads about 8 to 10 
years of data for ground water simulation. In absence of 
such a length of data, whatever data was available was used, 
Be;Core using the validated model for studying the anit- 
waterlogging measure sensitivity analysis on the parameters 
of the model was done* Four parameters namely S, T, Q and 
L were subjected to the sensitivity analysis. It was found 
that S,T are not very sensitive compare to 0 and L, It is 
therefore necessary to know the data on Q and L more care" 
fully than S and T. 

In all, mathematical modeling, determination of suitable 
boundary conditions have to be decided before validating and 
using the model. In the present case there was some doubt 
about the proper boundary conditions to be used. In order 
to varify the effect of different boundary conditions on the 
results, three types of boundary conditions v^ere studied. 
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These are (i) Constant head, (li) No flow and (iii) No flow 
on part of the boundary and specified flow on the remaining 
part of the boundary. It was found that these different 
boundary conditions do not have moh g^fect on the results 
when we move away from the boundary. How for away? can 
be determined by gradually reducing the mesh size. In the 
present case a luesh size of 6000 m x 6000 m was used. It 
was found that 3 k im. away from the boundary the affect of 
boundary conditions is negligible, Hov/ever one has to 
reduce the mesh size and recompube the results. This was 
not attempted in the present thesis for want of time. 

The validated model is then used for controlling 
the ground water table. It was decided to see how much 
numpage is needed to lower the ground water table 3 metres 
belov^ ground level wherever necessary. It was found that 
- 5 ^ 1 ^ ]yiCM of water per year is required to be pumped for this 
purpose. This water is of course available for irrigation 
purpose. It is the experience of the project authorities 
that during the summer months there is not enough water in 
the Kandana reservoir of the command area. Therefore, this 
extra water can be used for irrigation during suminer. Of 
course for meaningful results one has to apply the model not 
for a year as a whole but monthwrse or atleast for quarterly 
periods. This would require enormous amount of data and 
analysis. Due to lack of time this was not attempted. 
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A second policy ot lowering ground water table to 
6 m below ground level was also studied. This yields about 
MCM of water per year. This additional water can 
then be used not only in Summer but also for Rabi irrigation. 

b,2 SUGGESTIONS FOR FURTHER WORK 

Some of bhe limitations of the study have already 
been pointed in the above section. Further work can be 
directed towards removing these li.mitations. For example 
data should be collected atleast monthwise for 8 to 10 
years for ground water modeling. Mathematical modeling 
may then be done and validated using the split sampling 
method that is about half the available data is used for 
validdbing and the other half for checking. It may be 
necessary to develop two models one for monsoon period and 
another for non-monsoon period. The effect of boundary 
condition should be more thoroughly investigated using 
smaller grid intervals near the boundary, 

B'rom the present investigation it has been found 
that Q and L are more sensitive than S and T. It is therefore 
necessary to check and to see whether the available data on 
Q and L are quite reliable. Further the grid size used in 
the present model namely 6000 m x 6000 m is rather coarse. 

It may be necessary to reduce this to 1000 m x 1000 m. This 
would entail collection of more data but the results will be 
more realistic. 
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Finally it is necessary tc' update the model period!* 
cally by using the new data as and when they become available. 
It has been the experience of many ground v/ater modelers 
that the models requires changes as new data become 
available . 
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